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THEORETICAL ANALYSIS OF STATIONARY POTENTIAL 
FLOWS AND BOUNDARY LAYERS AT HIGH SPEED* 

By K. Oswatitoeh and K. Wieghardt 

The present report consists of two parts. The first part 
deals with the two-dimensional stationary flow In the presence 
of local supersonic zones. A numerical method of integration 
of the equation of gas dynamics is developed. Proceeding from 
solutions at great distance from the tody the flow pattern is 
calculated step by step. Accordingly the related body form is 
obtained at the end of the calculation. 

The second part treats the relationship between the dis- 
placement thickness of laminar and turbulent boundary layers 
and the pressure distribution at high speeds. The stability 
of the boundarj 1 " layer is investigated, 'resulting in basic 
differences in the behavior of subsonic and supersonic flows. 
Lastly, the decisive importance of the boundary layer for the 
pressure distribution, particularly for thin profiles, is demon' 
strated. 


PART I 
NOTATION 


p pressure 

p density 

T absolute temperature 

k ' ratio of specific heats 



*"Theoretische Untersuchungen dber stationare Potentlalstromungen 
und Guanos chich ten bei hohen Geschwindigkeiten. " Lilienthal- 
Gesellscha.ft fur Luf tfahrtforschung Bericht S 13 /l.Teil, pp, 7“2^» 
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H coefficient of friction 

w velocity vector 

w magnitude of velocity 

u, v velocity components 

$ velocity potential 

c velocity of sound 

c* critical velocity of sound 

Ma = w/c Mach number . 

F = 1 - Ma 2 

0 W ^ stream density 

Re* Reynolds number of the displacement thickness 

5 boundary -layer thickness 

f 

5* displacement thickness 

t 3 momentum thickness 

f stream filament section 

r radius of curvature of the stream line 

N normal to the stream line 

E maximum bump elevation 

R related radius of curvature 

Subscript n refers to the conditions in the free -stream region, 
subscript a to the outer flow, subscript w to the wall. Subscript 
m refers to the chamber or drum values in the phase quantities and 3n 
the velocity to the highest obtainable value. The quantities of state 
in part I are made dimensionless by the chamber quantities and the 
velocities by the highest velocity obtainable. 
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PART I 

1. ROTES ON THE CHARACTERISTICS OF COMPRESSIBLE POTENTIAL FLOWS 


The equation of gas dynamics is derived hy means of the energy 
equation rather than the adiabatic equation as customary. A simple 
formula is obtained for the stream density which is valid in a wide 
range about the critical velocity of sound. By applying this formula 
a simplified equation of gas dynamics is derived which in the tran- 
sition zone from subsonic to supersonic, for small velocity components 
v , describes the processes very accurately. Lastly, the problem of 
flow around a cylindrical body, symmetrical in two directions, is 
analyzed. It is found that, from a certain flow velocity on, located 
above the critical velocity, no maximum velocity can occur at the 
point of maximum body thickness. 

In the description of a gas flew the most general case involves 
six unknown functions, namely, the three components of the velocity 
and the three phase quantities of the gas, the pressure p, the 
density p, and the absolute temperature T. The equation of state 
of the gas permits the temperature to be expressed in terms of the 
pressure and density, thus leaving five unknown functions for the 
calculation of which the three Euler equations and the continuity 
equation are available. For the miss’ ng equation it is customary 
to use the adiabatic curve to eliminate the pressure and density from 
the equation ana so arrive at an equation between the velocity com- 
ponents, that is, the so-called equation of gas dynamics. However, 
it appears to be unknown that for the derivation of the equation of 
gas dynamics tho assumption of the adiabatic is not necessary at all, 
but that the U 30 of the energy theorem itself is sufficient. This 
derivation is briefly carried out in the following, while having 
recourse to the vector method. 

Pressure, density, and temperature are made nondimens ional by 
the corresponding "chamber quantities" p , p m , T m ; that is, the 
quantities of state at velocity 0, and all occurring velocities 
and velocity components by the maximum obtainable velocity, that 
is, the velocity at pressure 0. With w denoting the velocity 
vector, c the specific heat at constant pressure, and k the 
ratio of the- specific heats this maximum velocity is 
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The equation of energy of an ideal gas in stationary flow has 
expressed nondiraansionally - the following simple form 

*2 + |=l (1.1) 

and the corresponding continuity condition reads 

W f \ 

div v + -=■ grad p = 0 (1.2) 


The oft repeated quantity 

2 

n= 3T rT 

has a physical significance] it indicates the degrees of freedom 
of a molecule. For air n = 5 is very exact. 

The Euler eqtiation is then written as 

^grad — - w x rot ^n = - -~ grad p (1.3) 

The quantity xn enters the equation through the nondimens ional 
notation. 

Pressure and density are eliminated by forming the gradient of 
the energy equation, thus obtaining 

0 - grad + ~ grad p + (w 2 - 1) ~ grad p 

Scalar multiplication by w and application of (1.2) and (1.3) 
gives the equation of gas dynamics 


(l - v 2 ) div w + nw grad \ = 0 (1.^) 

This equation is written here in a form where the velocity of sound 
is already eliminated and only the flow velocity itself is present. 

Me will use this equation in the next section. 

For an insight into the potential flow properties in the speed 
range of Ma = 1 the just derived equation is much too complicated. 
So the processes in a flow filament are analyzed, unsteady variations 
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and friction processes excluded so that adiabatic changes of state 
can he assumed. The energy equation (l.l) then affords a connection 
between velocity and density, and the stream density vp can be • 
represented as a function of the velocity w: 

6( w) = pw = w(.l - w 2 ) 11 / 2 (1.5) 

It is icnovn that this function reaches a maximum at the point where 
the velocity is exactly equal to the velocity of sound. This 
particular point is generally denoted as the critical speed of 
sound c*. With f as the section of a flow filament the continuity 
equation reads fS = const. The speed c* 5s therefore characterized 
by the fact that a flow filament for this value of w reaches a 
smallest possible cross section. Tor v > c* as for w < c* the 
flow filament section is greater. 

Less familiar is the smallness of the stream density changes 

6 

over a very substantial speed range . To indicate it '£•* is 

represented for x - l-40‘ in the range of 0.5c* < w < 1.5c* in 
figure 1. Quantity denotes the va^ue of 0 for w = c*, the 

same applies to the derivatives of G- This characteristic of the 
stream density is of decisive significance for the effect of the 
boundary layer on the flow, as will be shown elsewhere. 

Consider the function 0 in the vjcinity of the maximum 
developed and signify its derivative with 0* w etc. Now 

it is found that the parabola 



is already sufficiently accurate for a wide speed range. This 
approximation is indicated by dashes in figure 1. The calculation 
for the coefficient of the quadratic term gives the simple result 


c* 

1 n *2 w k, + 1 _ n + 1 

e ~e*~~ ~ 2 n 


( 1 . 6 ) 


The equation (1.6) serves in good stead for the derivation of a 
simplified dynamic gas equation for two-dimensional flows on the 
limiting assumptions that the y component of velocity w, 
signified hy v, is small compared to the velocity of sound and that 
u, the x component of the velocity, does not differ too much 
from the velocity of sound. The stream density wp can be 
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replaced by up = ©(u) and the equation of continuity (1.2) on 
applying the same omissions as effected with respect to the terms 
with the factor v for the derivation of the Prandtl law, can be 
written as 


— Ji 

p ox 


(up) 


+ 


<Yv u 

by ~ Q 



+ 



= 0 


The coefficient of ^7 depends only- on u; it is simply a ^dif- 
ferent method of expressing the well-known quantity 1 - Ma - . 
Assuming, aside from the smallness of v, that jt- © u can .be 

regarded as constant results in the Prandtl-G-lauert .analogy. 

If this coefficient were plotted against u in the vicinity of 
the sonic velocity, it would show that it can assume negative as 
well as positive values and at u = c* is equal to zero. So 
the premise of constancy of this quantity can no longer be main- 
tained, especially since the derivative © u changes signs at 

sonic velocity, as seen from figure 1. The variations of H on • 

0 

the other hand are no more weighty than the variations of the 
entire coefficient anywhere in the range of not too high subsonic 
speeds. Thus In support of Prandtl* s law ~ can very well be 

put equal to this quantity in the free-stream region, but not for 
6 U . This equantity is computed by (1.6) and gives 


_ u o 8* r 
9 q c* 


+ 1 ) 



= 0 


The subscript o denotes the quantities in the free-stream region. 
Using the notation 


_ up e* , 

+ h (u - 

for 

1 and ~ ~ 
c*. 


- U 


(1 - 7) 

1 the simplified equation of gas dynamics 

6u dv 


rix by 


= 0 


(1-7) 


U can be positive and negative. Here also the introduction of 
a velocity potential is accompanied, although in simplified form, 
by the undersirable change of the equation from the elliptical 
to the hyperbolic type. To secure solutions which have supersonic 
zones by an analytical method It is advisable to find solutions of 
(l.7)j because it combines the simplifying assumption of small v 
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with a very accurate description of the processes in the critical 
sonic speed range. This was the reason for the brief derivation 
of the equation. 

The fact that the flow filament section has a minimum at the 
critical speed may, under certain circumstances, have very 
characteristic consequences for the velocity distribution at the 
appearance of supersonic zones on bodies, as will be demonstrated 
for the case of two-dimensional flow past a body that is symmetrical 
about two mutually perpendicular axes. The flow direction is to 
be along one body axis, that is, the angle of attack equal to zero. 

The flow is to be adiabatic and irrotational, the latter 
characteristic being expressed by 

dw w 
"cTn = ~ r 

where K is the normal to the streamline and r Its radius of 
curvature. The sign for N is so chosen that it is positive when 
the normal points out from the radiu3 of curvature. Equation (1.8) 
holds exactly for all two-dimensional potential flows. By the 
continuity condition in the form. 

Of = Constant 

and the freedom from rotation (1.8) the flow is completely defined. 
The origin of the coordinate system x rnd y is placed in the 
center of the body, axis x is made coincident with the flow 
direction (fig. 2), and tho aran oF positive, y value ana L.vzed . 

The cylindrical body is visualized as being exposed to a flow 
velocity which Heads to the formation of a supersonic zone near 
the thickest part of the body and it is assumed that in every 
stream filament the maximum velocity is reached at the point 
x = 0, an assumption which certainly should be fulfilled for 
subsonic flows. A point on the y — axis with supersonic speed 
must have a maximum stream filament width, a point with subsonic 
speed, a minimum of stream filament width. In the supersonic 
region the curvature of the streamlines on the positive portion 
of the y— axis must decrease less rapidly than on concentric 
circles, in the subsonic zone the stream line curvature must 
decrease more rapidly- than for concentric circles. Hence no. 
great error is introduced when in the vicinity of the point on 
the y— axis where sonic velocity is reached, the streamlines are 
replaced by concentric circles, and it will not lead us far astray 
when this is assumed up to a value of y equal to twice the 
distance from the cylinder of the point with the critical sonic 
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Telocity c*. After the streamline curvatures are approximately 
known the velocity distribution on the y— axis in this zone is 
completely defined by (1.8). If the piece which the body cutB off 
from the y— axis is denoted by H and the radius of curvature of 
the profile on the y_axis by R, its velocity distribution is 

u R r-, rs\ 

= R — H + y (1 * 9) 

According to (1.6) it may be stated that the volume of flow through 
a section of the y_ axis is then greater than on an identical 
section of the free stream, if at a _,.a.r ti< ula? ^ oint the inequality 


u Q < u < 2c* - u Q 

is fulfilled. Bince the velocity distribution for x = 0 is 
defined by (1.9) up to the constant u g., it also is the 
difference in through, -flow volume for "7 > H in the free— stream 
region and on the y— axis. It may now be asked at what value of the 
constants the absolute amount of this through flow difference reaches 
its highest possible value and the answer is found in the fairly 
accurate equation 


Uy= H = 2c * ~ u 0 (l-9a) 

that is, that the stream density on the y— axis must nowhere be 
less than in the free stream. For a simple picture it is imagined 
that (1.9) with the constant (l.9a) is applicable up to the attain- 
ment of speed u Q and that from this y value on, the constant 
flow velocity prevails. This break may occur at the value y = y Q 
for which the equation reads 

As near the body more can flow by than on a strip of equal width 
in the free-stream region, because of the increase in density , we 
must proceed from the cylinder only as far as the free stream is 
displaced. The result is therefore a highest possible value of H, 
denoted by which is given by the equation 

9 o-^max = J (® “ ®o^ (l .10) 

y ~^max 
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The integrand is given "by (1.6), (1.9), (1.9a), H is to he replaced 
by since u is equal to w on the y — axis. The evaluation 

of (1.10) gives the following relation between flow velocity and 

®max 

E 


TABLE I 


^0 

c* 

0.70 

0.75 

0.80 

O.85. 

0.90 

0.95 

1.00 

H 

R 

0.053 

0.026 

0.013 

0.0059 

0.0020 

0.0004 

0 . 


The extension of the speed by pieces at y = y Q unquestionably 
introduces an error; but it can only cause a. shift in II - while 

not changing the existence of such a value. In the subsonic zone a 
streamline may be regarded 'as a bump and the residual rise in 
through flow volume due to increase in velocity computed by an approxi- 
mation process that applies In the subsonic range. The result then 
is a finite variation of the integral in (1.10) and a correspondingly 
different H^^. The possibility of a speed increase in y direction 
in the subsonic region must be rejected, as it would invalidate the 
present considerations. Hence it is seen that the assumption cannot 
be applied to all bodies and therefore the following principle: 

To each flow velocity u Q , there corresponds a definite ratio 
H ffl^ /B. If the ratio H/R exceeds this limit for a body symmet- 
rical in two mutually perpendicular directions and lying along the 
flow direction, there is no flow for which velocity maximums can be 
reached cn the entire y — axis. 

It must be expected that the maximum speeds on the y — axis 
disappear only in the supersonic range . But since It cannot be 
assumed that velocity maximums in the supersonic range disappear 
on a part of the y—axis while a velocity maximum appears on the 
body, we are led to the following principle . 

From a definite value of H/R on, for bodies and flow directions 
of the described type, there is no flow at which a speed maximum with 
local supersonic zone Is reached at the point of maximum thickness 
of the body. 

A boundary point for these specific values of H/R is given 
In table 1. 

In the subsonic zone this principle has no analogy. 
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2. METHOD FOE THE NUMERICAL INTEGRATION OF THE EQUATION 

OF GAS DYNAMICS 

A numerical - graphical method is indicated for finding 
solutions of the equation of gas dynamics with supersonic zones, 
hy progressive calculation of the entire flow, starting from an 
exact solution at great distances from the body. The exact body 
form follows at the end from the shape of the streamlines . Exceeding 
the sonic velocity causes no special difficulties or peculairities. 


Limited to two-dimensional, irrotational flows with w = grad $ 
past a cylindrical body, equation (1.4) gives for the velocity 
potential cp a nonlinear differential equation of the second order 


D 



1 


(n + 1) 



n 

dx 2 


+ 


1 



(n + 1) 



& 

2 

ay 


( 2 - 1 ) 


2n 


0$ 8$ 8 2 $ 


The zero point of the coordinate system is placed in the body, its 
dimensions are of the order of magnitude of unity, and the flow 
strikes the body along the positive x-axis . The boundary conditions 
Tor $ then read: 




= 0 at the body itself, N 


§-*> 


and 


denoting the normal, and at infinity 
— >u = (dimensionless) flow 


for z = yx 2 + y 2 —>c 
velocity. 

On passing through the local velocity of sound 



equation (2.1) changes from the elliptic to the hyperbolical type. 

For this case the exact integration has been successfully secured 
for single specific examples only. For the subsonic range several 
general approximate solutions are available, the simplest and best 
known of which is the solution <J> obtained by the Prandtl rule. 

This satisfies equation (2.1) better as the body becomes more slender 
and the distance from the body becomes greater. 
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The following method is therefore indicated. Compute the 
’Prandtl solution 0_ for the entire flow and attempt to secure 
the correction cp in such a way that 0 + 9 = <±> becomes a 

solution of the complete equation (2.1). As the analytical 
calculation of cp is too complicated, a numerical method is 
advisable, starting from the outside (z » l) where cp = 0, 
and progressively continuing inwardly toward the body. The 
exact body shape follows at the end of the calculation from 
the streamline distribution; however it is to be suspected that 
it essentially remains similar to the form of the Prandtl solution. 

For this purpose the differential equation for 9 = 0 — 0 
is set up; is an exact solution of the considerably simplified 

equation (2.1;: 


with 


(1 - Ma 0 2 ) 



+ 




( 2 . 2 ) 


The subscript r - rofers to the conditiona at infinity, cp 

fulfills the complete equation (2.1) up to an error which 

can he computed by means of (2.2): 



- 2n 


Ox 3y 3x 8y 


(2.3) 


Putting 0 = Op + 9 in (2.1) and regarding it as a differential 
equation for 9, ' Cp follows as term of zero degree in 9. As 
9« Op in the entire range, it is assumed that it applies to the 
derivatives as well. Merely the terms of the zero and first degree 
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and the greatest term of the second degree in cp need to he included. 
There results 


V W _ 3y 2 


" £p ~ L VSJy “ W J a? 


+ 21 Ma„ + n 


8 2 0p a&p _ 8 2 $ p 8^ 89 


ox' - 8 x 


dx by 8 y cix 


j n , o 1 8 2 $>p c-3)p c 

+ 2 J(n 1) Ma 0 2 ~ n J §rT zf + n “ 


8 2 $p 8 § p [ 89 

dx 8 y 8 x i Sy 


8 % oOp 8-9 _ , . . 8!>p ^9 89 

8y ox dy 8x ^ x 2 ox 


(2.4) 


The tetter 1 0 p satisfies the equation of gas dynamics, the 

easier is the determination of 9 . Go, at great distances from the 
tody the equation can he substantially simplified. For z » 1, 

O^n 8'9n 

especially on slender todies, : hence we can put 

oy ox 

= 0 and — = 0 in 2.4 hut in contrast to the Prandtl rule 


consider the variation of 


Since 9 0 for z->co, the term 


of the second degree is omitted also. And equation (2.4) is simplified 
to 


14 = - fa - Fn + F 


'!)■ 


% ,a . 


(2.5) 


1 - (n 4. 1 ) f 


si- Ma^ = 


84uV 
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and 


F» 




2 


This equation could also "be derived from (1.7); F is an abbre- 
viation for — 0 U . 

The boundary condition for z 00 is cp = 0 where the Prandtl 
rule applies exactly. On the other hand, however, the disturbance 
of the flow by the body is very extended when the flow velocity 
approaches sonic velocity. It is therefore necessary to determine 
an initial approximation for cp analytically so as not to be com- 
pelled to start at unduly great distancee. For this purpose (2.5) 
is transformed further. While Ma = Ma 0 , by the Prandtl rule, hence 

/ d&p \ 

F = F_, the more exact term F = Fh + F’ . ( u ) , so that 

o o V 8x 


F + = _p» 

0 8x2 Sy 2 


8 2 cp\ f ^V 

vox 2 Sx 2 / \>, x °) 8x 2 c)x 


S 2 ®P 8cp 


(2.7) 


There 


8 2 cp 

8x2 


8 2 ( 


can be ignored with respect to 




8x 2 


and it is assumed 


8cp 8&-> 

that « ■=— - — u 0 at great distances; for example, in the cal— 
ox 8x u 

culated case: For a parallel flow and dipole, — u Q dies down as 

l/z. 2 , but ^ as l/z^. Hence finally: 

8x 


+ t, 

o. p v 2 v 0 ’ 8*2 

o~£r oy ox 


( 2 . 10 ) 


All the equations for the process are now available. The general 
process of calculation is as follows: First determine the entire 

field of flow of the incompressible fluid numerically, and then the 
compressible approximate solution by the Prandtl rule for a fixed 
Ma number of the flow velocity. The correction cp on a strip for 
great y follows f?om (2. 10). From here on 9 is computed 
numerically, step by step. At great distance from the body we 
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therefore first use equation ( 2 . 5 ) and later in the neighborhood 
of the body the more complete equation (2.4). Having thus 
determined cp for the entire field, the stream lines and there- 
fore the body contour itself, as well as pressure and velocity 
distribution, are obtained from § = <$p + cp. Naturally the process 
can be built on another approximate solution; however, the formulas 
probably become simplest when the Prandtl rule is used. 

For the present the range of application of this method is 
confined to the flows where the Prandtl rule affords a good 
approximate solution and the fundamental assumption cp << $ is 
actually fulfilled. 1 


Excluded are accordingly flows around not sufficiently slender 
bodies, as well as areas in which the velocity of sound is sub- 
stantially exceeded. This also manifested itself in the calcu- 
lated example. Flow around a cylinder (circular in the incompressible 
case) at Ma = 0.7454. The calculation was perfectly smooth into 
the supersonic range, where it had to be broken off especially 


since ^5* an & 
ox 


3 2 C> 


8 2 cp 
3x ciy 


quickly rose to the order of magnitude of 


£ 

3x 


and 


3x 3y 
difficulties. 


but exceeding the sonic velocity itself involved no 


In principle we can also free ourselves from the approximation 
that cp « $p, when in the formulation of equation 2.4 we consider 
terms of the third degree in cp; the length of the calculation, 
however, becomes disproportionately large. In another more appro- 
priate method the assumption cp « is omitted and the tedious 


calculation of §p in the entire field of flow is eliminated. The 
previously described cp method is utilized only for computing the 
initial values for large Z. The new method is as follows: fu is 

evaluated at large distances from the body for y » y^; where yq 

is chosen so large that the error of the prandtl solution is suf- 
ficiently small; cp is evaluated from equation 2.10. This affords 

the exact solution of the dynamic gas equation $ = 1 $^ + cp in an 

initial strip. From here on § itself is calculated step by step. 
The width of the initial strip from —xq to + xg 1- must extend 
upstream and downstream from the body so that for all y at 
x < — xq and x > xq the Pranatl rule is applicable with sufficient 
accuracy. From y-, on, where as is then known, 

a 2 $ 

— - — yq 


3y 2 


is graphically extrapolated to 




for certain 


■“•For bodies which are symmetrical relative to the y-axis 
also, Xq is naturally = Xg. 
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fixed abscissae x (Ay is the length of the stepj it can he 
assumed quite large at first, and reduced again later in proximity 
of the hody) . Next 


8$ ! 


_ S$ _ S 2 $ 

" 8y2 


Ay 


yr A y yi “ yi-^y 


is plotted against x and graphically differentiated, which gives 


S 2 $ 

5y~5x 


^2$ 

Plotting t .- 2 . , - against y and integrating gives 
dy ox ' ~ ' 


y "Ay 


s$ 

Sx 


S$ 

Sx 


y x -Ay 


y l" Ay s5 

n 8 c o 


rsE 37 


yi yi 


Lastly the variation of 


SO 

s 


over x yields 


differentiation. With it 


yi 


S g o 

Sx 2 


hy graphical 


SO SO S 2 0 S£& 

S? Sy ^ STSP Sx 2 

are known for y = y - Ay. 

From the equation of gas dynamics (f-’.l), in which the simpli- 
fication ~ = 0 can he made so long as it is valid that 
oy 

is calculated for the required x values and plotted 

Sx ap 2 

against y. Then the calculation is repeated, extrapolated 

for yq “ gAy and so forth. 

S 2 0 

If the values of — — computed for y > y^, the value 

Sy 

extrapolated at y, - ±Ay and that computed for y, - Ay do not 
j. 2 j. 

form a smooth curve, the step must he repeated with a differently 

extrapolated value for the particular abscissae . In this event it 

is hatter to reduce the length of the step. Since the differentiation 


i 
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of the curves 


8 $ 

ST, 


and 


is uncertain at the 


y=const . 


y=const , 


■boundary points x = — x-^ and x = x 2 , it is advisable to compute 


$ = + <p also in two vertical strips 


x < - x-j_ and 


x > x 2 


and 


to Join the progressively defined points to these edg^ strips. 


The direction of integration for this step method must be 
chosen at right angles to the flow for the following reason. At 

flow around a body exposed to a flow along x. — — is sure to be 
8(5 ox 

greater than r— almost everywhere, especially in the supersonic 
oy aa& 

zone. Thus at entry in a supersonic zone the coefficient of — - 

8x2 

in (2.1) goes through zero. This does rot interfere in the above 
method since ( 2 . 1 ) is used for computing 

If, however, we integrate in the x ““direction and solve 
equation ( 2 . 1 ) for 
coefficient of 


82$ ° x 


8y c 


then difficulties will result. The 
can, on the other hand, disappear only far 
8 $ 


above the speed of sound where — 3.3 not important at the point 

dy 

under consideration. We can also prove this state of affairs with 
the help of characteristics. When discontinuities in the velocity 
or their derivatives appear we cannot integrate across a character- 
istic. On the other hand the characteristics of our flow become 
nearly vertical so that again we can not calculate in the x direction 
in this supersonic region. 


3. ILLUSTRATIVE EXAMPLE 


A flow symmetrical in x and y is computed for a Mach number 
of flow of 0.7^54. The flows on smooth bumps with supersonic zones 
are obtained exactly, but on the other hand the flow past a closed 
body is obtained only with errors in the region of the stagnation 
point . 


The described method is tried out on a very simple example; 
we start from the incompressible flow (subscript i) past a circular 
cylinder 


= 


Xj + 


*i 


2 P 

*1 + n d 


(3.D 
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the free-stream velocity and the radiua are taken as unity. Prandtl's 
rule is applied to a fixed Mach number of flow Ma c to which the 
dimensionless flow velocity 


3$ l 
u o = 33T/ 


Ma, 


Ma + n 


Z — ^±00 

corresponds. The abscissas for this transformation are contracted 
by \jl - Ma 0 2 . 

(3-2) 


X s X 


fl - Ma Q x ± 


P “ \l 2 

The ordinate^ remain the same: y ~ y^ - so that 


3l> 




Sr ' “o SJ 


0)2^ u 

p _ o 


d'di. 


S/x - Ma' 4 ^ 


(3-3) 


To compute cp by equation (2.10) the coordinates Xj and y = y^ 
are used, so that 

2\ -^cp S 2 cp _ 3 2 cp ^2^ 


(l - Ma 2 ) -Ll , Z2 - $1$ 
V ° / ^v2 + 


3x d 3y d dx-^ 3y t_ 


u 


1 - Ma, 


- 1 + Ma ^ 

o 


3x^2 


1 - (n + 1) u Q 2 (dj^/dx^ 1 


1 - u r 



(3.*) 


Development of the right-hand side for large ts \jx±^ + y~ by 

using equation (3*1) gives for the first approximation, when 
. 2 ? 


£ = 


x. 


n = “2 = 1 


and 
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£p(xi, y) 


! *Lq!!.. — Sg. (1 - ft, + 8,2) 

nu 0 \Jl - Ma 0 2 z i 



( 3 . 5 ) 


A particular integral of this Poisson equation is obtained 

with the help of the separation formula cp = — pf (£) . Thus the 
general solution is written zq 


cp = 


— 12 Ma r 


nu 0 \/l - Ma Q S 




+ *pot 


( 3 . 6 ) 


with Acpp ot = 0 and c arbitrary. 


As boundary condition for cp 
it shall be small compared to £> p , 

than ~r ~ . But for the rest 


Vt 


the sole requirement is that 
that is, decrease more rapidly 

and c can be chosen at random. 


The physical meaning of this ambiguity is as follows: Owing to the 

disregarded terms of higher order in the formulation of (3-5) only 
the effects of the first order of the body at great distances are 
taken into account. But these are the same for different body 
forms. So the calculation yields different section forms, depending 
upon the choice of c and Tpot- The manner in which c and 
<p ^ affect the body form cannot be evaluated until several examples 
hive been worked out . Up to now only one such example has been 
worked out, owing to lack of time. 


The Mach number of flow Ma 0 = |/ 5/9 = O.jh-^b had been 
specifically chosen. The dimensionless velocity is then 

u G = \jl/ 10 = 0 . 3162 and x = x p = — - Xq . 

In (3.6) only half cp p 5 0 G = 1/8 were assumed for 

simplicity, thus eliminating the linear term in £ 
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cp = 




(3.7) 


with this <p the velocities (and their derivatives) of the exact 
solution $ = + 9 for y = 10 and 0 < x< 6, as well as to 

x = 6 , 0 < y< 10 were computed and $ was determined for 
y < 10 and x < 6 by the described step— by— step method. In view 
of the symmetry of flow relative to x and y the calculation 
in one quadrant was sufficient. The step length Ay up to y = 1.5 
was Ay = 0.5; from there on 0.2. 


While the exceeding of the sonic velocity (first at y = 1.35) 
caused no difficulties, the calculation could not be carried out to 
the body because of another reason but had to be broken at y = 0.6. 
For at x K 0.6 the horizontal components of the velocity 


(39 

changes so rapidly for smaller ordinates y that the graphical 

differentiations became too uncertain to compute the next step. 

As is seen from the contour of constant velocity (fig. 3) a 
steady but still very sudden rarefaction occurs and on a point 
symmetrically situated with reference to x = 0 a compression 
occurs. This phenomenon would of course not be plain at a lower 
flow velocity, but it is certainly characteristic of the flow in 
proximity of the stagnation point where the speed increases quickly 
from subsonic to supersonic. For this point of the flow field 
another method must therefore be developed. 


So while unable to obtain the flow around a finite bodj r 
with a stagnation point, the data obtained thus far are nevertheless 
very informative for subsonic flows with supersonic zones. The 
calculated streamlines and lines of constant velocity are shown in 
figure 3. Visualizing, the lowest streamline in figure 3 as rigid 
wall, we get the flow along a smooth bump with a supersonic zone 
near the highest point. Since this streamline is already very 
steep for x^ 0.65* it can be assumed that the velocity distribution 
around the finite body (with symmetry axes x = 0 and y = 0.5) 
indicated in figure 5 is fairly accurately reproduced by the dotted 
line. Incidentally, it is rioted that even Prandtl’s rule yields 
considerable errors near the stagnation point . 

Figure 4 shows several streamlines magnified five times in 
elevation, along with the respective velocity distributions. Not- 
withstanding the similarity of the individual peaks the velocities 
differ considerably at various places. The velocity — and with 
it the pressure distribution of thin bodies — is therefore at high 
flow velocities markedly dependent upon the exact shape of the body. 
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Noteworthy also is the steep velocity increase at a point where 
the streamlines themselves are still comparatively flat. 

The contours of equal velocity in the supersonic region prove 
the principle set up in section 1 according to which the highest 
speeds under certain assumptions do not occur at the point of 
maximum thickness of the body. Even the equation (1.9) applied 
for the derivation of this principle is satisfactorily confirmed 
in figure 6 , where the velocities on the y— axis are plotted along 
with the hyperbola (dotted) that touches the curve w(x = 0 )/c* at 
w/c* = 1. From the far-reaching agreement of the curves it follows 

that in the vicinity of w = c* the expression u = w(x = 0 ) <= — - — 

b + y 

is a good approximation. With this example the accuracy of table I 
can be checked. In view of the flow velocity of u c /c* = 0.7746 
H mp /R = O.OI 9 would have to be expected according to this table, 
but by the calculated example it is proved that from H„ n y /R = 0.031 
on, the speed maximum is no longer situated at the greatest ordinate. 
Thus, it is seen that table I is a good representation of the order 
of magnitude of H y /R. The difference is attributable to the 
fact that the hyperbola used for the approximation gives too low 
speeds in the subsonic range. 


P A R T II 

4. INTRODUCTORY NOTES ON BOUNDARY LAYERS AT HIGH SPEEDS 


Studies of the behavior of supersonic flows in parallel 
channels disclose that In the supersonic zone, principal flow 
and boundary flow are in unstable equilibrium In certain circum- 
stances. An effect of the boundary layer on the principal flow 
in the zone of the critical speed is to be expected for the reason 
that here small variations in stream density cause considerable 
changes In speed. This is particularly plain in the calculation 
of the flow through a Laval nozzle at high subsonic speed with 
observance of the boundary layer. 

In order to gain an insight into the condition of the 
boundary-layer flow at high speeds, which we will study in the 
following, consider an example from the sphere of Incompressible 
flows, where the conditions are better controlled. We consider 
the circulation— free, incompressible, and stationary flow around 
a circular cylinder at a high but still subcritical Reynolds 
number. Computing the pressure distribution at the body with the 
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aid of the potential theory on the assumption that the cylinder has 
no dead -air region 'behind it and then calculating on the basis of 
this the boundary-layer conditions, say, with the aid of a refined 
Pohlhausen method, we find a separation point in the zone of rising 
pressure. It is found that the omission of the dead-air region was 
wrong. The pressure distribution on the body must therefore be 
computed with due allowance for the dead-air region and then it can 
be hoped to attain a result corresponding to reality when the 
dead-air region is so assumed that the related pressure distribution 
yields separation exactly at the starting point of the free stream- 
line. This example shows that potential flow and boundary -layer 
flow usually depend upon each other. In general, we can say that 
.the potential flow determines the boundary -layer flow, also that 
the boundary -layer flow determines the potential flow. The former 
can be stated with great approximation in flow without pressure rise. 

It is a known experimental fact that for large expansions a- ■ 
flow simply does not follow the boundaries of the regionj but it 
should be remembered that for the development of a dead -air region 
not the expansion of the stream filament but the fact of a lores sure 
rise is decisive, which only in subsonic flows goes hand in hand 
with an increase in stream filament section. In supersonic flow 
on the other hand a contraction of the stream filament results in a 
pressure rise. Thus visualizing a parallel channel with a flow of 
Ma > 1 a too strong growth in boundary layer caused by some 
disturbance is followed by a pressure rise, which in turn favors 
a stronger growth in boundary layer. In contrast to subsonic flow, 
an unstable equilibrium of boundary layer and principal flow is 
involved in this instance and a very considerable boundary layer 
growth must be reckoned with in certain circumstances. It may, 
in a straight channel result in a sudden strong pressure rise at 
the flat wall and so in the formation of a dead-air region (fig. 7 (a)). 
(Compare reference 11.) If the pressure rise is so great that the 
flow becomes subsonic, the relation of main -flow and boundary -lay er 
flow is stable again, the dead-air space cannot remain in this part 
of the channel. If a small pressure rise is involved of, say, a 
weaker oblique compressibility shock, the principal flow experiences 
a directional change in the sense of a channel contraction. The 
dead-air space must increase wedge-like, but this holds only over 
a short distance, otherwise the flow would have to revert into the 
subsonic range. It is therefore' to be assumed that at an oblique 
compressibility shock, as met with in figure 7 (a), the turbulent 
Intermingling imposes a limit on the growth of the dead -air space. 

These qualitative reflections lead to the conclusion that in the 
range Ma > 1 an unstable state of equilibrium, must be reckoned 
with in certain circumstances between principal - and boundary -layer 
flow, which may promote the formation of dead -air regions even at 
a flat wall. 
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A disturbance of the unstable stato of equilibrium of principal 
and boundary-layer flow in the supersonic range is favored by the 
fact that any minor disturbance in a supersonic flow is propagated 
undamped along Mach lines. Thus, a pressure rise in a supersonic 
tunnel can be dispersed by a small disturbance far upstreamj on the 
other hand, the pressure rise sets in again some distance downstream 
as we can also infer from our example. 

The unstable behavior of the boundary layer in the supersonic 
zone must disappear when the principal flow approaches sonic 
velocity. In the critical speed range w = c*, which is of particu- 
lar interest in flows past bodies with high speed, the fact stands 
out that this is the range of maximum flow density. But the proce- 
dure in computing the incompressible flow past an airfoil is such 
that the pressure distribution is obtained from the potential, flow 
without consideration of the displacement effect of the boundary layer, 
and then the boundary layer is computed with the aid of this pressure 
distribution. This is not permissible however in the region of the 
velocity of sound, because a minor variation in stream density 6 
exerts a very substantial effect on the speed. This is readily 
apparent in figure 4, where peaks with comparatively minor form 
changes produce very unlike pressure distributions. This of feet 
increases with increasing flow velocity* 

The effect of the boundary layer on the flow in the vicinity of 
the velocity of sound is illustrated by a simple example, which, 
although it involves no flow problem, is novertheless informative 
for the appraisal of the displacement effect of boundary layers 
at high speeds. The velocity distribution in the nozzle used by 
Stanton (reference l) for his experiments was computed by appli- 
cation of the simple flew filament theory, once without boundary 
layer, and once on the assumption of 'a laminar boundary layer. 

The boundary-layer calculation is made with the help of a process 
which will be explained in the following section. The initial 
value of momentum and displacement thickness at x = - 0.20 was 
estimated. The dimensions of the nozzle are so small that it can 
be assumed that no turbulent transition takes place. Stanton’s 
test series C is illustrated in figure 7(b). The velocities were 
determined by measuring the static pressure on the axis of the 
axially symmetrical nozzle (lower test points) and adjacent to the 
wall (upper test points). The theoretical curves by Oswatitsch and 
Rothstein (reference 2) and the flow filament solution with and 
without boundary layer allowed for are included for comparison. 

The former was computed only as far as the separation point. It is 
seen that the asymmethy ’is reproduced qualitatively correct by the 
flow filament solution with boundary layer taken into account. 

The displacement thickness at the narrowest point of the nozzle is 
not quite 2 percent of the nozzle radius. Computing the velocity 
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distribution for the same nozzle in incompressible flow with and 
without consideration of the boundary layer, the results in both 
cases are essentially even lines. Even at speeds about 15 percent 
below those of test C, any boundary-layer effect is quite insignificant. 
This may be taken as proof that the asymmetry in nozzle flows which 
at the most, manifest local supersonic zones, are caused by boundary- 
layer effect . As to making the computation, only the following is 
mentioned. That one gets at first the distribution of the displace- 
ment thickness from the stream filament solution and then a new 
stream filament solution taking into account the calculated displace- 
ment thickness is proof in itself that such an iterative procedure 
is permissable at very high subsonic speeds. Displacement thickness 
and stream filament solution are obtained step-wise at the same time 
in the downstream direction. 

The influence of the boundary layer on a submerged body will be 
handled in section 7. Our example, however, shows that we cannot 
hope to obtain results that correspond to the real process in some 
degree, for the flow problem with high velocity, without examining 
the boundary layer. We then have to remove, in practice or experiment, 
the influence of the boundary layer, perhaps by suction. 


5. CALCULATION Off DISPLACEMENT THICKNESS OF LAMINAE AND TUEBULENT 

COMPEESSIBLE BOUNDARY PAYEES 


For more accurate calculations on boundary-layer effect in 
flows at high speed, formulas for the variation of the displacement 
thickness are necessary. This is done in this section, first for the 
laminar and then for the turbulent boundary layers . In the derivation 
of the formulas for laminar boundary layer a refined Pohlhausen method 
is given for computing laminar compressible boundary layers for given 
pressure distribution. 

Inasmuch as the pressures in the compressible zone transverse 
to the flow direction in the boundary layer is also to be regarded 
as constant, the relative density variations within the boundary 
layer are in amount equal to the relative temperature variations. 

Thus, if sonic velocity prevails in the outer flow the relative density 
variation inside the boundary layer amounts to about 20 percent, 
since stagnation point temperature can be approximately assumed at the 
wall. So, at not too high supersonic speeds a qualitatively identical 
behavior in the boundary layer and in the incompressible range Is to 
be expected. 


24 


NACA IM No. 1189 


So as not to exclude the possibility of compressibility shocks 
beforehand, the flow outside of the boundary layer is called outer— 
or— principal flow instead of potential flow. 

Being primarily interested in the beliavior of the displacement 
thickness, the behavior of other quantities is studied only to the 
extent that it appears in the result without loss of time. This is 
the case in the study of laminar boundary layers, where the momentum 
thickness is comparatively easily obtained and the displacement thick- 
ness derived from it. 

The process is based upon an improved Pohlhausen method in 
conjunction with the reports by Bohlen (reference 3 ) and Walz 
(reference 4) . 

Boundary-layer equation and continuity condition for the 
stationary case are written as follows 


cSu du dp b / 3u\ 

pu 35E + pv ^ = ~ d3E + Sy 3yJ 

( pu) + !y ( pv) = 0 


(5.1) 


The coordinate system is chosen in the usual manner so that x is 
tangential and y normal to the contour of the body. Thus, in the 
general case, x, y imply no Cartesian coordinates in the following: 
p is the friction coefficient dependent on temperature. The 
quantities p, p, u, and v are not made dimensionless. It further 
is assumed that the boundary— layer thickness is small relative to 
the radius of curvature of the wall, so that curvature effects can 
be disregarded. The second Navier-Stokes equation gives then 
exactly as for incompressible flow the result that p is merely 
dependent on x, but not on y, which in the boundary— layer 
equation was already evidenced by formation of the ordinary 
derivative of p. 


After integration of the boundary-layer equation over y, the 
application of the continuity condition gives the von Karman momentum 
equation which is written in the form 


d_ 

dx 




du a 
p a u a dx 



(5.2) 
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The subscript w indicates the Talues at the wall, subscript a 
those in the outer flow. In the momentum equation displacement 
thickness and momentum thickness are defined by 



(5.3) 


where 5 is the so-called boundary— layer thickness, which is now 
chosen so great that 5* and can be regarded as independent of 
8. The displacement thickness has the physical significance that 
the through flow volume in the boundary layer is reduced by an 
amount that, on the assumption of pure potential flow, is equivalent 
to a shift of the wall by piece 8* in positive y direction. 


Equation (5-2) can also be given the form 


u A 2ai! + ad! 

a dx dx 




(5.2a) 


Aside from the usual assumptions with the aid of which the boundary- 
layer equations are derived, no restrictions of any kind have been 
made so far. Only the problem without heat transfer on the wall, 
the so-called thermometer problem ^ is tree ted in the following. 
Further we will make one approximation, by which we will specify 
the form of the boundary— layer profile by only one parameter in 
addition to the Mach number of the outside flow. Next it is 
necessary to make an assumption concerning the configuration of the 
density profile. Having seen that in the vicinity of the critical 
velocity the velocity profile especially might be decisive, while 
the density variation is unimportant , the case of laminar boundary 
layers is limited to the assumption that the temperature at the 
wall always attains the tank temperature T m and satisfies the 
energy theorem within the boundary layer. This ties in .also the 
assumption that temperature — and velocity — boundary layers are 
of equal thickness. Accurate calculations on flat plates indicate 
that this assumption also holds in a considerable supersonic range 
(reference 5) • 

The pressure in the y direction being constant, the density 
variation follows from the temperature variation as 
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The parameter for the boundary— layer characterization is derived from 
the known houndary conditions which for u = v = 0 in the boundary- 
layer equation loads to 


cp. 


w 


Su 

oy 


w 


+ - 


du 

p a u a dx 


Since the internal friction is only dependent on T and no heat 


transfer takes place, 


[ j = 0, hence with application of & 

V Sy A 


X* = - 




.HI) 


u 

Uc 


w 


p a ^ 


du„ 


P w dx 


(5-5) 


We have avoided the introduction of the boundary— layer thickness 
itself in the equations other than at the unimportant place as upper 
limit of the integral. The version of (5-5) was largely taken from 
Walz’s report (reference h) . The parameter X* differs from the 

conventional Pohlhausen parameter by the factor (— V, the density 

\ & / 

refers to that outside, the internal friction to that at the wall. 
For a class of velocity profiles, such as for M = 0, for instance, 
the individual profiles which are represented by the magnitude of 

parameter X* the quantities and f ~ ~ \ 




and I — — \ can be taken as 

\ & / y 

from the class of profiles, hence 


function of the parameter 

p 

q P a & 

obtain u a by (5-2a) simply as function of 

0 X U-y 


X*. Choosing 


the Pohlhausen profiles as profile class gives the curve of Bohlen— 
Holstein (reference 3) in figure 1, while the Hartree or Howarth pro- 
files result in Walz’s curves of figure 1. Moreover, it is not 
necessary at all to have an analytical representation of the profile 
class, it can equally be given as family of experimental curves. 

If the outer flow is dependent on a Mach number, one profile class 
is used for each Mach number. In order not to come to grief because 
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of our ignorance in the sphere of compressible velocity profiles 
two known facts are takon advantage of: First, we know that the 

velocity profiles on the plate at constant pressure are not very 
closely related with the Mach numbers of the outer flow (reference 5) • 
Therefore, this is assumed to be the case in retarded or accelerated 
flow also; secondly, we know that the single parametric method in 
the incompressible, which utilizes the Pohlhausen profiles, leads to 
fairly practical results, although the Pohlhausen profiles themselves 
do not represent the actual profiles very well. An exact represen- 
tation of the velocity profiles themselves is not needed, the main 
point is the displacement thickness for which the integration over 
the profile form is already accomplished. 

On the basis cf these arguments we therefore select, independent 
of the Mach number, the Hartree profiles for the velocity distribution, 
which, as a rule, probably represent the incompressible boundary 
layer, best of all. The density variation is then given by (5*4) as 
function of the velocity distribution and owing to the presence of 
u a 

quantity — as function of the Mach number. 
w m 

ft-# 

By a mimerical method the q-uantities and 

V 

then obtained as functions of A* and Ma, as exemplified for the 

u~ 

following Ma which correspond to the — values : 


Ma 

r — i 

0 : 

1 : 

1.2 

1.5 

2.0 

^a 

w m 

, 

! 

0 1 
i 

_L 

o.bQo ! 
1 

0.473 

' 0.557 

O .667 


We obtain the curve system of figure o. 

The curves are shown at the left as far as the separation 
point, at the right they proceed to the point up to which the 
Hartree profiles are calculated. Several curves were extrapolated 
beyond it, and indicated by da.shes . The curve Ma =0 is identical 
with the Hartree curve from Walz . 



Thus with the velocity distribution, the variation of the Mach 

2 


number of the outer flow, and an initial value of 


a 




the uarameter 


A* can be formed; with it and observing Ma the quantity 
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u 


can be taken from figure 8 and the variation of * 

l J w 


a dx |_i 
"w 

computed. From it we obtain again this quantity at a point shifted 
by one step and the calculation can then be repeated (reference 3 ) • 

Q nd 

6 * - « ■? + i.tq cs vi-via , Pov»vi/*^ 4- r* nw-m -i +■ /T\ ^ 


Although 


is wanted, it was preferred to compute 


Mv 


because the equation for this quantity is very much simpler. Know 
ing the outer flow, S ^ can now be specified . To determine 5* 

C u 

thus further requires ~ which is a function of A* and Ma, 


represented in figure 9 ■ 


Since the quantity in the boundary -layer equation is to 

dx 

be defined as accurately as possible the following formula is of 
advantage 


_1_ dS* = j3_ _o_. 5* dA* ±_ _o_ /B*\ JL_ % 
8* dx 6* dA* ■Q dx 5* u a \d / dx w m 

s — - 


+ 1 / 1 jL p ^ 2 + Ha£ 

2 [p a -a2 dx ^ ua 
\ M-w 



(5-6) 


This formula contains only quantities dependent either on Ma and 
A*, or which can be taken from the previous boundary -layer calcu- 
lation. The derivative with respect to — . which is a function 

w m 

of Ma, was preferred over that with repect to Ma for reasons of 
simplicity. The first term at the right-hand side in (5*8) is 
generally the principle term. 


If it is desired to eliminate 


dA* 


d 


, and , 
dx dx 


P ^_2 


in (5-8) 


so as to secure merely in relation to — — , -- -- - (which 


dx 


dx' 


dx ,: 


jt, j/, 

then because of ~~ enters in the equation ) and coefficients 
dx 


solely dependent on A*, Ma, and Re, the expression becomes 
fairly long, but since this relationship is used later it is 
given here, the Re number being suitably referred to the out 
side density, the coefficient of friction at the wall and the 
displacement thickness: 


Re* 


u aPa s * 

hw 


(5-7) 
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The formula reads 


as* 

Ox 


1 /s*\ 2 / _S jA 6* 

Ee* Vd J u u a to 


p /fL u \ j, 5 S* 

\aj ">j» 8 * * 


6 * 

•S 


2 + Ma d 


U& ±_ _a_ S* 
+ w m s * aHa 13 


+ 





&*£ d2 % 

u a 2 a*2 


5 * 2 

u a / 



+ 4 



a l _ S* du a 
Re* 2 u a to 


a-^Re* 




to 2 



(3-8) 


It contains the first and second derivative of u a made dimension-- 

less with the displacement thickness and the outside velocity , and 
also Re* and coefficients that are dependent on A* and Ma. 

The coefficients cc.-^, a?, And 0,3 for A* = 0 are given in the 
following table: 


TABLE II 

LAMINAR BOUNDARY LAYER } A* = 0 


Ma 

! 0 
1 

1 

x 1 

1.2 

1-5 

2-0 

^a 

i 

! 0 

! 

0.408 i 

0.473 

0-557 

O .667 

w m 

j 

1 



a l 

j 1-47' 

2.02 j 

2.30 

2.76 

3.82 

a 2 

: 5 . 0 I 

l 

4.89 1 

j 

4 .So 

3-87 

2.37 


1 0.3^9 

0.284 1 

1 

0.256 

] 

0.219 

0.154 

— 


Later on the equation is to he applied to the case -where the 
velocity distribution differs little and monotonically from 


Ua = Const, so that 


5* £ua 


B*f d2 % 
^ to 2 


are also regarded as 


and 
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"being small . Restricted to the linear terms in the derivatives of 
u a , the term with can "be struck out and the coefficients ctg 

and ctq taken at the point X* - 0 for the specified Mach number. 
The dependence of ap on the derivatives of u a proves to be so 
small that the same can be done for this quantity too . 

So on the assumption of small derivatives of u a simply 
(5-8) is taken with the constants of table II for the corresponding 
dS* 

Ma, for — of a laminar boundary layer. 

CL.X- 

To obtain a formula for the variation of the displacement 
thickness of a turbulent boundary layer a different procedure is 
required. On analyzing the cause of the variation of an incom- 
pressible turbulent velocity profile at a specific pressure 
variation it is found that the pressure forces are primarily 
responsible . The shearing stresses introduced by the turbulent 
intermingling play, however, a subordinate part. It is true that 
the difference of the two effects is not so far reaching that a 
second profile could be computed accurately enough from the 
specified velocity profile when the shearing stresses are discounted, 
because the shearing stresses are able to substantially modify the 
character of the profile ; but for the calculation of the variation 
of displacement thickness, which essentially involves an integral 
over the velocity variation in the incompressible case, the shearing 
stresses can be disregarded. 


Hie result at Ma = 0 is the following approximation formula 
for the turbulent boundary layer: 


dS* 

dx 



(5.9) 


As the integrand is always positive, it can be taken from this 
formula that a speed increase is accompanied by a decrease in 
displacement thickness and a speed decrease by an increase in 
displacement thickness. At constant outside velocity the displace- 
ment thickness remains constant, according to (5*9). This result 
is naturally wrong, as indicated by experiments on the plate at 
constant pressure. For in this case the variation in displacement 
thickness is contingent upon the turbulent shearing stresses, so no 
correct result is to be expected. The formula could be improved 
by the addition of the conventional formula for the variation of 
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the displacement thickness, hut it would serve no useful purpose, as 
will he seen. It is of greater significance that in contrast to 
equation (5-3) for the laminar boundary layer the second derivative 
of u a is lacking in (5.9)* Since (5-9) was obtained by several 
rougher omissions its practicability is illustrated in figure 10. 

The experimental values of — and are 3hown plotted 

u a dx dx 

against the arc length x of Gruschwit' 1 ’, ’ s (reference 6) test 
series 3, along with the variation in displacement thickness calcu- 
■ / v de- 
lated by (5 '9), the integral being formed at = 0. It is found 

that the formula reproduces the actual conditions adequately, as far 
as the area of greater accelerations, where errors begin to be intro- 
duced. This is, of course, due to the fact that 5* =0 imposes a 

limit on the decrease in displacement thickness . 

These experiences in the incompressible zone can now be 
interpreted to the effect that the turbulent shearing stresses for 

the calculation of - can also be cancelled in the compressible 

zone. But even this assumption is insufficient to develop a law 
for the variation in displacement thickness; additional data on 
the density distribution in the boundary layer are needed. In the 
case of turbulent boundary layers the energy theorem i3 not 
directly applicable, because the density -boundary layer Is 
probably twice as great a3 the velocity -boundary layer; hence, the 
density varies in an area in which the velocity is already practically 
considered constant (fig. ll) . The result of it is that the varia- 
tion in density plays the same role in the calculation of 6* as 
the variation in speed within the boundary layer. Unfortunately 
only one measurement of a turbulent supersonic profile is available, 
and naturally there is little sense in developing a theory without 
further basiB. However, in order to reach a tolerably correct 
numerical value, the part of the boundary layer in which the density 
alone varies is disregarded for the present, since it involves 
only about 10 percent of the displacement thickness, and, in the 
remaining portion, putting the stream density as a function of the 
velocity as follows: 


pu 

p a% 



(5-10) 


H to be taken from the oxperiment. Now the derivatives of p can 
be expressed by derivatives of u, 1^, and p a with the aid of 
(5*10). This enables us to derive a formula for the variation of 
displacement thickness, neglecting the turbulent shearing stresses. 
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as* 

dx 


5 * iu§, 1 
% 






(5.11) 


H' is the derivative of the function H according to the argument 

— and S the -place where -!L can be put equal to unity 
% ~ % 

= lj , while S* represents the correct displacement thickness 

hence, integrate up to a point where itself is equal to unity 

Pa, 

/-E- = l\ . This means, we state that the 10 -percent of the dis- 
\Pa ' 

placement thickness between the point — = 1 • and = 1 

% Pa 

contributes to the variation of the displacement thickness an amount 
which corresponds to its portion of the displacement thickness. 

For Ma = 0, equation (5»ll) naturally changes to (5-9) • If the 
density and speed in the boundary layer are specified, the integral 

can be evaluated also . We have calculated the expression in paren- 

du a 

theses for a profile by Gruschvitz, for which =0, and for the 

velocity and density profile represented in figure 11; thus we obtain 
the constant for two values of the Mach number. 

TABLE III 

TURBULENT BOUNDARY LAYER 

! Ma | ~0 ; ITT) 

1 ; ,1 

I <x 2 | 5.1 2.2 | 

The close agreement of coefficient for the turbulent and the 

laminar velocity profile is noteworthy. 

6. STABILITY STUDY ON THE FLAT PLATE 


A study of the equilibrium of boundary layer and supersonic flew 
on the flat plate indicates that an -unstable state is involved. The 
growth of a small disturbance in a laminar boundary layer differs 
somewhat from that in a turbulent layer and is, especially in the last 
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case, very, rapid. In incompressible flow a stable equilibrium 
exists between principal flow and boundary layer. 

Having secured the variation of displacement thickness 5* 
in relation to the velocity variation of the outer flow, the 
reciprocal effect of principal ~ and boundary layer flow is now 
analyzed in the simplest case, namely, in the flow at the plate, 
without specified pressure distribution. 

Since the effect of Bma.ll disturbances is to be involved, the 
Mach number of the outer flow Ma is regarded as constant and. the 
v component of the velocity considered small relative to the 
velocity of sound. After introduction of a velocity potential the 
simpld equation (-2.2) is involved, and written in the form 

' ' Oyy = ^ (Ma 2 - 1) (6.1) 


The x-axis is to be in plate direction, the y— axis normal to it. 


Now it is necessary to represent the effect of the boundary 
layer on the potential flow in form of a boundary condition. The 
boundary layer is therefore visualized as being replaced by an 
elastic layer superimposed on the plate, which has the property 
of always attaining the thickness equivalent to the displacement 
thickness of the boundary layer at the particular place for the 
prevailing velocity distribution. That is, the equation 


v 


u m 
dx 


must satisfied for y = 5*. 


(6.2) 


This condition is inconvenient to the extent that the boundary 
for which it is to be fulfilled is not specified beforehand. But, 
inasmuch as the disturbances are to be small, hence the outerflow 
is to differ very little from a flow u a = Const., the boundary 
condition for displacement thickness &* in undisturbed flow is 
assumed. By assumption the departure of 5* from the value of 
the displacement thickness for the undisturbed flow must be small. 
Hence it seems immaterial whether v is specified at y = &* or 
at y = &* + db* in the linearized problem. Besides, the study is 

.. as* 

to be restricted to such a small area that itself can be 

regarded as constant at u a = Const. 
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Furthermore the boundary condition (6.2) has the property of 
giving the same v component of the velocity at y = 5* as a 

dB* 

boundary layer with equal on the assumption of potential 

flow in the entire space . 

In (6.2), v and u are none other than the components of 
the outer velocity, hence in the notation of the proceeding section 
equal to v a said u a . 

a e-#. 

Applying (5*8) or (5. 11) to — — - gives then as boundary 

condition of the problem, linearized in the derivatives of u, the 
following equation for 


y = 6*: v = u _ a 5* 8u _ a r 0 * J^u s * 2 (6.2a) 

Re* 2 J 5 x 2 

If a laminar flow is involved bhe corresponding constants 
must be taken from table II } if, turbulent flow, table III; in the 
latter case, must be put =0. In view of the linearization 

B* and Re* must also be regarded as constant although the 
variation of u in the first term is not Important, it is considered 
nevertheless, because the solution then is reduced to the treatment 
of a homogeneous linear differential equation, which means some 
simplification. 

Now it is attempted to find the solution for the case that 
the plate Is exposed to a flow with the velocity u(x,y) = u Q = Const 
and at a point x = 0 at the plate the velocity is artificially varied 
by an amount <$C u Q . The coordinate system is turned through a 
small angle so that the x~axls in point x = 0 is exactly in flow 
direction and the y«-axis normal to it. The tangent of the angle 
of rotation is defined by the variation of the displacement thickness 
at x = y = 0, which is equal to 


dB * _ _ ( ^1 
dx " Re* 

Strictly speaking a transformation of the coordinates in the 
equations themselves should be effected. But since the boundary 
layer itself makes no difference between these two directions, 
and so a rotation merely involves more paper work without any 
physical significance, it is disregarded and the equations applied 
to the new coordinates. The coordinates are in addition visualized 
as being made dimensionless by the displacement thickness and the 
origin shifted to the point x = 0, y = 8*. 
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These new coordinates are denoted with 


x* 




and after introducing the velocity potential in (6.2a) give the 
following differential equation with the respective boundary 
condition 


9 y*y l = ®x , x' (Ma 2 — 1) (6.1a) 

at 

y* = 0; = - a 2 ®x*x' ~ ^^x t x , x t (6. 2b) 

Assuming a very general solution of (6.1a), and writing the 
potential as sum of a potential of a principal flow u Q and a 
small disturbance 


$= u 0 8*x‘ + u&* f^x* — \j Ma^ — 1 y’^) 

+ g^x* + \Jmh 2 - 1 y’j | 

f and g are arbitrary functions of which it is merely required 
that their sura at x 1 = y* = 0 be equal to unity. It is seen 
that g gives Mo.ch lines which point toward the boundary layer, 
hence stem from a disturbance from the outside. This function is 
thus put identically zero since such disturbances are to be 
disregarded. Introdpction of the thus obtained solution in the 
boundary condition gives the functional form of f . Denoting the 
derivative with respect to the argument 


T) = x* — Ma 2 — 1 y* 

with subscript q, we get 

$x* = u Q S* + $x»x* » l»*f . • 

^x’x’x 1 = ub*^^; $y f = - u8*^Ma 2 — If 

which inserted in (6.2b) gives an ordinary differential equation of 
the form 

a.Ke*f + a 0 f — y Ma 2 — If =0 
3 Wl 2 Tin v tj 


(6.3) 
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This equation is easily solved. First postulating a laminar boundary 
layer, hence ^ 0. Then from the requirement for x' = y 1 = 0: 
f = 1, a requirement for the second derivative f can be satisfied, 

because the upper equation can be regarded as differential equation 
of the second order of f p . Since at this point only the consequences 

of small velocity disturbances, not the consequences of disturbances 
of the velocity difference are to be studied, the added requirement 
for x' = y* = 0 is f = 0, which gives the solution 


f 


n 




eVl 


where t^ and t 2 are abbreviations for the expressions 


a. 


“ 1,2 


111 ~2 \ , VMj 

■ W 4 \ aoSe*/ a 


'' 2 


VMa 2 - 1 


2 a^Re* ^ W 4- ^a^Re*/ a^Re* 


(6.4) 


As Re* in general has the order of magnitude of 1C)3, the last 
term under the root is a term of greatest influence . Thus for 
appraisals at high Re* we can put 


t 


1,2 


♦v 

A; 


h Ma 2 - 1 
1 a^Re* 


(S .4a) 


It is to be noted however that the critical Re* which corresponds 
to a value of about 1.4 x 10^ must not be exceeded as will be 
shown later. 

By use of (6.4) the velocity distribution in a laminar boundary 
layer on the plate is obtained as : 


u 


5* 


v “ Ur . 


U 


t l " t 2 


t 1 e t 2 


^x 1 -^Ma 2 - ly - t 2 e^l (x'-\jui 


■5a 2 -! y'J 


(6.5) 


which by (8.4a) is reduced to the simple form 


f \/Ma*~ ~ 1 


u ft u Q + u oosh\j^~- (x 1 - \/Ma2 - 1 y j 


(6. 5a) 
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If a turbulent boundary layer were involved, hence a, = 0, 
the first summand in ( 6 . 3 ) cancels out and only one boundary condition 
can be satisfied. Again requiring f = 1 for x* = y* = 0 gives 


f 


n 



The velocity field under the assumption of a turbulent boundary layer 
at the plate is 


u = u D + ue 



( 6 . 6 ) 


From ( 6 . 5 ) and ( 6 , 6 ) It is seen that the boundary disturbance along 
Mach lines is propagated into the flow. The Interference velocity 
u is always accompanied by a function which grows considerably 
with rising value of the argument, while in the case of the laminar 
boundary layers the coefficient c, plays the principal part. In 
turbulent boundary layers the coefficient a g is essentially 
involved. Thus the boundary layer of a flat plate in flow with 
constant velocity is in both instances In an unstable state of 
equilibrium with the principal flow, which with observance of the 
terms of the first order only, lets a small disturbance grow 
infinitely. The type of growth is, of course, quite dissimilar 
on the two boundary layers. To secure a measure for the Instability 

of the state, we may ask for which value of x f =-g-#r at y* = 0 

the disturbance has grown to twice the amount and call this quantity 
the length of growth A. It is not made dimensionless by the 
displacement thickness. 

The length of growth in a laminar - boundary layer A^ Is assessed 
by (6.5a), The hyperbolic cosine grows for a value of the argument 
of around 1.3 to the amount 2. Accordingly 


A, - 1.3. 


a-^Re* 


i 


Ma^ — 1 


5* 


(6.7) 


The length of growth of a turbulent boundary layer Is 


At 


0.70 




( 6 . 8 ) 
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Postulating a laminar boundary layer at Re* = 1000, tables II 
and III give the, following length of growth 


TABLE IV 


Ma 

1.2 

1-5 

1.7 

2.0 

A 7 

8* 

25 

18 


12 

At 

8* 



1.1 



Noteworthy is the unusually short length of growth in the turbulent 
boundary layer; but even that in the laminar layer is still very 
small when bearing in mind that the displacement thickness in 
supersonic flows is of the order of magnitude of 10“3 to 
10“ 2 centimeters. 

The investigation was restricted to small disturbances. The 
extent of growth once they have reached greater amounts remains to 
be proved. One thing is certain that the outerflow cannot increase 
to great velocities, because the boundary layer cannot drop below 
the amount 8* = 0. Thus no limit in velocity decrease is imposed. 

It may be presumed that the velocity decreases until the boundary 
layer breaks away. In general, the instability of the discussed 
equilibrium condition will become evident in a pressure rise, 
probably an oblique compressibility shock. It would not be sur- 
prising if oblique compressibility shock occurred in the center on 
a flat wall (fig. 7(a)), The example cited here could be multiplied 
by many others, perhaps even by flow around conical tips. It should 
be kept in mind that a pressure rise can cause transition of the 
boundary layer. In the example adduced here the boundary layer is 
already certainly turbulent. 

This study of plate flow can be regarded as first result in 
this sphere of instability of Supersonic boundary layers. It would 
be desirable to get away from the assumption of small disturbances 
and constant flow velocity. This seems altogether possible by a 
combination of characteristics method and boundary layer computation. 
For the turbulent boundary layer, of course, the laws of variation 
5* would have to be analyzed first. 

One unusual fact is that in the measured pressure distribution 
on a wing, such as those by Gothert (reference 7), for instance, 
pressure increases were almost never observed in the supersonic 
zone, except in form of compressibility shock or occasionally at 
small Reynolds numbers, where laminar boundary layers must be assumed. 
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It appears entirely possible that this fact might he explainable 
by the cited properties of the supersonic boundary layer. 

The corresponding behavior of a laminar boundary layer in 
incompressible flow (Ma = 0) is briefly indicated. The disturbance 
at great distances from the wall, that is, for great values of y, 
must disappear. On these premises, (6.1a), (6.2b) by the same method 
of calculation give 

u = u D + ue cos (p 2 xl “ 017*) (6.9) 


with the abbreviations 



The decisive term at high Re* numbers is again ct^Re*. For 
Re* = 500 it approximately is 


Pi 


P 2 


1 

^ 2cLyRe* 


0.053 


that is, a strongly damped oscillation is involved. The analyzed 
equilibrium of laminar boundary layer and outer flow in the sub- 
sonic zone is extremely stable according to it. This method of 
analyzing offers the further possibility of exploring the stability 
of laminar subsonic boundary layer relative to nonstationary dis- 
turbances and comparing the results with Tollmien t s calculations 
(reference 8). For nonstationary velocity variations Poh.lhausen r s 
method is, of course, not practical in general, in the form given 
here . 


Incidentally, the requirement of damping of the disturbance 
for great y is not fulf illable in subsonic flow on the assumption 
of a turbulent boundary layer at the plate . This result may have 
its cause in the fact that (5-9) does not meet all requirements. 
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7. SIGNIFICANCE OF BOUNDARY LAYER IN THE PRESSURE 
DISTRIBUTION ON A BODY 


Appraisals indicate that the flow in the critical range of 
sonic velocity is very substantially affected by the boundary 
layer. Without its inclusion a correct calculation of the pressure 
distribution therefore seems, in general, not very promising. In 
many instances the behavior of the boundary layer actually governs 
the pressure distribution. 

On examining the pressure distribution at a bump computed in 
section 3, (fig. 4), a symmetrical velocity distribution is also 
found on a body symmetrical about the y-axis. This is, however, 
in great contrast to the experience in tests (compare, fig; 12) , 
where symmetrical peaks were invariably accompanied by asymmetrical 
velocity distributions . Naturally the guest ion is whether there 
is only one solution for each bump but it will be shown that, owing 
to the boundary— layer effect, symmetrical solutions can be expected 
as little as in the example of the velocity distribution in a nozzle 
(fig. 7(b)). 

By (5*8) the displacement thickness of a laminar boundary layer 
for constant outer speed is 


5 * 



p w 

u a p a x 


What is the possible extent of the bump in order that the boundary 
layer remain laminar? Figuring with tests in a low-pressure tunnel, 
the values at critical velocity are 


u a = 3 X lO^cm/sec ; p a = 0.8 X 10"3g/cm3; p w = 1,8 X 10 _J+ 

CGSE 


It is to be presumed that the critical Reynolds number at sonic 
velocity does not differ substantially from that in incompressible 
flow. Taking the critical Reynolds number formed with the plate 
length at 


Ee crit . 


= 5 x 10 5 
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gives the critical Reynolds number (5.7) formed with the displacement 
thickness at 


Re* 


crit . 


1.4 X 10 3 


with the previous values of Ug., p a , and the critical values of 

plate length and displacement thickness are . 

V 

x crit . = 3-8 cm; &* crit . - 1.1 x 1CT 2 cm 


So in order to prevent transition from laminar to turbulent flow in 
the boundary— layer model., lengths of only a few centimeters may be 
permitted in the usual test arrangements, provided that no strong 
accelerations are involved. 


Conversely, the critical length indicated here gives a measure 
for when the transition point is to be expected on a plate flow in 
an exhaustion tunnel at sonic velocity. In a free-air test this 
length is reduced by about half because of the higher density. 

In the schlieren photograph of an infantry shell in flight at 
around sonic velocity (fig. 13 ) (references 9 and 10) the oblique 
compressibility shock is evidently released by transition, its 
effect being probably amplified by the unstable behavior of the 
boundary layer. The fact that a missile at small supersonic speed 
is involved is immaterial; since a straight compressibility shock 
prevails in front of nose of the missile, it actually flies as if 
in a subsonic flow. 

Analyzing the bump in figure 4, which at the point of its 
greatest height has nearly constant sonic velocity for some distance, 
and supposing the points of strong velocity rise and velocity 
decrease (x = 0.6) to be about 2 centimeters apart, the displace- 
ment thickness at the peak is certainly greater than that of a plate 
1 centimeter in length in flow at sonic velocity. Therefore 


5 *x=0 > °*56 x IQ" 2 cm 


At the point of substantial speed decrease, separation must be 
definitely expected. A calculation by the expanded Pohlhausen 
method shows that the momentum thickness grows with increasing arc 
length. Much greater is the rise in the ratio of displacement 
thickness to momentum thickness (fig. 9) which for Ma = 1 increases 
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from point X* = 0 to the separation point from value 3.2 to 4.7. 
Considering the fact that the momentum thickness itself increases 
up to the separation point, an empirical rule can be established 
according to which the displacement thickness is doubled between 
X* = 0 and the point of separation. 

The difference between the displacement thickness at the 
separation point and at the highest point of the peak is in the 
example; therefore 


6* ^ - 5* ~ O.56 x 10-2 cm 

separ. x=o 

The difference in height of the highest point and at point of 
separation 4 se p ar is (compare fig. 4) 


h ~ 3' x 10~2 =m 

separ. 

While the variation in h„„ ov . duo to the boundary— layer effect 
amounts to a mere 20 percent, the illustration shows that a change 
in height of bump by this amount must be followed by an extra- 
ordinarily great change in velocity distribution, so that there 
can be no question of attaining symmetrical results in the 
experiment . 

The conditions in the presence of a turbulent boundary layer 
ere considerably worse. A little calculation on Grushwitz’e test 
series 3 (reference 6) discloses that the displacement thickness 
multiplies from the point of transition to the point of turbulent 
separation by about 25 times. ^Assuming turbulent separation at 
the point of severe velocity drop the greatest displacement effect 
(height of bump + displacement thickness) would also exist on a 
bump of considerably greater absolute dimensions at the point of 
separation due to boundary— layer growth. It is supposed that the 
displacement effect of the body, increased by the displacement 
effect of the boundary layer, undergoes no substantial increase 
behind the highest point of the bump. In turbulent boundary layer 
and thin profiles or low bumps this is possible only to the extent 
that a compressibility shock occurs at the point of greatly reduced 
profile thickness; furthermore, a compressibility shock would have 
to occur so much farther downstream as the bump or the profile is 
flatter. It also is feasible that the effect of the increase in 
displacement thickness is raised by strong return flow behind the 
point of separation. These qualitative results can be checked 

against the work of Gcthert (reference ( ) . 
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The fact that a compressibility shock can occur when there is 
enough space available for the increaoed displacement thickness 
caused by it is to be regarded as reason for the fact that the 
separation computed by stream filament theory in figure 7(b) is almost 
exactly coincident with the start of the compressibility shock in the 
test. 


It may be asked how the streamline pattern in a flow problem 
must look, in order that the compressibility shock be possible. This 
can be answered to the effect that the compressibility shock on slender 
bodies is to be expected near the point of vanishing streamline curva- 
ture. Since the streamlines in the zone of critical sonic velocity 
are approximately parallel, the points of vanishing streamline curvature 
must lie near a common orthogonal trajectory, hence, a potential line. 
Along it the velocity changes little according to (l;8). In a flow 
that differs little from the critical sonic velocity, the free -stream 
velocity is therefore to be expected in the vicinity of points with 
zero streamline curvature. If the curve decreases rapidly at a place 
with supersonic velociby a decrease to the outer velocity must be 
counted on. The marked velocity variations in figure k coincide with 
the streamline inflection points. On flat profiles a point of separa- 
tion can be regarded as starting point of a free streamline with very 
little curvature . The streamline curvature must thus decrease very 
substantially in the separation point and it is seen that a strong 
compressibility 3hock produces through the separation connected with 
it a streamline pattern that favors the appearance of the compressibility 
shock. This argument is therefore not suitable for finding the location 
of a compressibility shock. 


8. CONCLUDING REMARKS 


The preceeding work shows that in a calculation conforming to 
reality the pressure distribution of a body in a flow at supercritical 
free— stream velocity may not be given by the potential flow, that the 
boundary layer plays a decisive role here. In general, the potential 
flow around the body permits not even an approximate calculation of 
the boundary layer. This means that in contrast to incompressible flow 
the pressure distribution on flat bodies can also be much different. 

It is therefore intended to first improve the process of calcula- 
tion of the potential flow with a supersonic region. With the process 
we will ascertain the flow around a substitute bod.y. This will have 
approximately the same displacement effect that is found on an 
experimentally investigated body including its deadwater region and 
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the displacement effect of its boundary layer. We can also anticipate 
from our calculation a strong velocity increase at the body nose and 
a strong velocity decrease at the point where the curvature of the 
substitute body disappears. 


Translated by J. Yanier 
National Advisory Committee 
for Aeronautics 
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Figure 1.- Relationship between stream density ' 0 /e * and velocity w/c* 

with approximation parabola for w£c*. 
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Streamline with the asymptote y = 0.6 
Streamline with the asymptote y = 0.8 
Streamline with the asymptote y = 1.0 
Streamline with the asymptote y = 1.2 
Streamline with the asymptote y = 1.5 






Figure 5.-- Velocity distributions along a body 
with local supersonic zone (extrapolated near 
stagnation point). 
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Figure 7b.- Velocity distribution in an axially summetrical nozzle according 

to Stanton (reference 1); test c. 

• measured on the axis 

o measured near wall 


— stream filament theory without boundary layer 

stream filament theory with boundary layer 

according to Oswatitsch and Rothstein (reference 2) 
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Figure 11.- Turbulent compressible boundary layer: velocity and density 
profile measured on the flat plate at Ma = 1.7. 
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Figure 13.- 



Infantry bullet at slight supersonic speed (reference 9 ). 
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